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SUMMARY 


The  purpose  of  this  paper  Is  to  present  some  approximate  methods  for 
finding  the  characteristic  values  of  the  integral  equation 

Xf(x)  -  f'  .»  coa  (x-*>  *  b  coa  f(y)dy 

0 

in  the  cases  where  a  or  b  are  small,  or  where  both  are  large. 
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A  NOTE  UN  A  CLASS  OF  INTEGRAL  EQUATIONS 
RELATED  TO  THE  BESSEL  AND  MATHIEU  FUNCTIONS 

Richard  Bellman 

31.  Introduction. 

In  this  note  we  discuss  the  behavior  of  the  solutions  of  the  integral 
equation 

(1)  MU)  -  J"  „a  c05  co»  (*♦/>  rtjrjdy 

0 

for  either  a  or  b  small,  and  for  both  a  and  b  large. 

If  a  ■  b,  we  obtain  the  equation 

(2)  Af(x)  -  ^  e2a  *COS  x  c°3  y  f(y)dy 

0 

which  we  recognize  as  Whittaker's  integral  equation  for  the  even  Mathieu 
functions,  [2],  while  a  -  -b  -  ik  yields  an  integral  equation  satisfied  by 
both  the  even  and  the  odd  Mathieu  functions. 

If  either  a  or  b  is  equal  to  zero,  we  obtain  an  equation  whose  charac¬ 
teristic  values  are  the  Bessel  functions  of  order  n  with  the  trigonometric 
functions  as  characteristic  functions. 

For  general  real  values  of  a  and  b,  the  equation  is  one  which  occurs  in 
a  problem  arising  in  the  statistical  mechanics  of  molecular  chains,  cf.  Montroll, 
[lj.  In  this  problem,  as  in  many  other  "nearest  neighbor*  problems,  only  the 
largest  characteristic  value  is  of  interest. 

In  the  second  section  we  use  the  simple  solutions  for  the  case  b  -  0  to 
obtain  power  series  developments  of  the  characteristic  functions  and  values 
for  small  b.  The  coefficients  in  the  power  series  developments  are  obtainable 
step-by-step.  For  the  case  where  a  and  b  are  both  small,  this  procedure  yields 
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results  which  seem  superior  to  those  obtained  by  merely  expanding 

(3)  ,a  c°a  (*-y)*b  cos  (x*y)  .  i  .  [a  c08  (x.,)  ,  „  cos  (x*y)]  *  •••. 


In  particular,  for  the  case  of  Mathieu  functions,  this  method  would  seen  to 
yield  some  new  expansions. 

In  the  third  section,  we  discuss  a  heurintic  method  for  obtaining  the 
asymptotic  behavior  of  the  largest  characteristic  value, 

§2.  b  small.  * 

It  seems  difficult  to  determine  the  radius  of  convergence  of  the  power  series 
we  determine  despite  the  fact  that  it  is  cle;  ~  a  priori  that  for  each  a  there 
exists  a  non-zero  radius  of  convergence.  Consequently,  we  shall  assume  that  b 
is  in  each  case  sufficiently  small  to  ensure  the  convergence  of  the  series. 

We  set 

(1)  \  •  XQ  ♦  \jb  ♦  X2b2  ♦  ••• 

f(x)  -  fQ(x)  ♦  fx(x)b  ♦  f2(x)b2  ♦  •••  . 

Substituting  these  expansions  in  (1)  of  §1  we  obtain  the  system  of  integral 
equations 

U)  XQf0  -  /  .*  cos  <*-  *)  f0(y)dy 

.  .0 

Vo  4  Vi  m  °q  **  008  (x’y)[V^  cos  4  fi(y)]<iy» 


and,  generally, 


(3) 


k*/-N  1 


2n 


/ 


a  cos  (x-y) 


Z  fv 


k*/-N 


cos  (x*  y) 
/! 


dy,  N-  -  0,  1,  • 


The  characteristic  functions  of  the  first  integral  equation  are 


U)  4<0).‘  2«I0(a), 

4m)  ■  2ni«(a)' 
« 


40)  -  l/VS 

f&}  -  sin  mx/Vn  , 
-  cos  mx/Vr  , 


where  IB(a)  is  the  Bessel  function  of  imaginary  argument. 

Let  in  what  follows  Xq  and  f^  represent  the  m-th  characteristic  value  and 
function  respectively.  To  determine  f^  we  use  the  fact  that  the  integral  equation 

1 

(5)  Mx)  -  J*  k(x,  y)a(y)ay  ♦  b(x)  , 

0 


where  X  is  a  characteristic  value  of  k(x,  y)  has  a  solution  if  and  only  if  b(x) 
is  orthogonal  to  the  characteristic  functions  associated  with.X.  Since  the 
equation  for  f^  is  as  given  in  (2),  we  obtain  for  X^  the  simple  expression 


(6) 


fQ  ®a  COS  coa  (x^  y)f0(x)fQ(y)dx  dy 


2n  - 

/  fo(*)dx 

0 


> 


which  may  be  evaluated  explicitly  in  terms  of  the  functions  I  (a). 

m 

Having  determined  X^,'  we  nay  solve  for  f^  by  means  of  its  Fourier  expansion. 
We  then  proceed  to  determine  X^  and  so  on. 

The  method  is  certainly  one  involving  a  good  deal  of  labor  if  higher  terms 
are  required.  However,  at- each  step  only  one  equation  in  one  unknown  need  be 
solved.  For  the  case  when  a  ■  b  and  b  is  small,  it  does  seem  superior  to  the 
process  depending  upon  the  equation  of  (3)  of  §1. 


§3.  Asymptotic  behavior. 

As  mentioned  above,  in  the  problems  encountered  in  the  mathematical  treatment 
of  physical  systems  involving  "nearest  neighbor"  interactions,  only  the  largest 
characteristic  value  is  important. 

The  method  we  present  below  is  frankly  heuristic.'  We  hope,  however,  that 

in  furnishing  a  guide  to  the  answer,  it  may  furnish  a  guideto  the  solution. 

Let  •  ■  aT,  b  ■  0T,  a,  3  real  and  positive,  and  let  T  »  1,  We  have 

'P  T(a  cos  (x-y)  ♦  3  cos  (x*y)) 

-  J  ®  f(y)dy . 

0 

Here  \  is  the  largest  characteristic  value  and  f(x)  is  the  associated 
characteristic  function.  Since  the  kernel  is  positive,  we  know  that  A  is  positive 
and  that  f(x)  is  positive. 

Setting  x  -  0,  we  obtain 

(2)  Kt( 0)  -  /n.T(a*^co.  y  f(y)dy< 

0 

♦ 

To  evaluate  the  asymptotic  behavior  of  the  right  side  as  T— >oo,  we  use  Laplace's 
technique  which  consists  of  observing  that  as  T — >-oo,  the  behavior  of  the 
integral  will  be  governed  more  and  more  by  the  behavior  of  the  integrand  in  the 
neighborhood  of  y  ■  0,  since  cos  y  has  a  unique  maximum  at  this  point. 

This  argument  is  not  rigorous  since  f(y)  also  depends  upon  T,  and  we  cannot 
be  sure  that  this  dependence  may  not  upset  the  previous  reasoning.  Having  issued 
our  caveat,  we  proceed  heuristically  and  obtain 
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> 


2,  T(a»U)[l  -  X-  .  ...] 

(3)  *f(0)  -  J'  e  0(0)  *  •••]  <*y 

f(a^)  2n  -y2^ 

■  f(0)e  P  e  dy  ♦  ••• 

0 

T(a^)  ^  -y2T (a*3)/2 

'■  f(0)e  y  e  dy 

0 


T(a*3)  ^ 
~  f(0)e _ yn 

V^T(a  ♦ 


Hence,  as  T  — ^od  , 

T(a^) 


(4) 


tp* T(a  ♦ 


provided  the  above  analysis  is  legit imate. 

For  the  case  where  3  ■  0,  the  above  is  the  well-known  asymptotic  expression 
for  I0(T). 
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